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Abstract 

The study of neutrinoless double beta decays of nuclei and hyperons require the calculation of 
hadronic matrix elements of local four-quark operators that change the total charge by two units 
AQ = 2 . Using a low energy effective Lagrangian that induces these transitions, we compute these 
hadronic matrix elements in the framework of the MIT bag model. As an illustrative example we 
evaluate the amplitude and transition rate of S~ — > pe~e~, a decay process that violates lepton 
number by two units (AL = 2). The relevant matrix element is evaluated without assuming the 
usual factorization approximation of the four-quark operators and the results obtained in both 
approaches are compared. 
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I. INTRODUCTION 



Lepton number violating (LNV) interactions with AL = 2 are widely viewed as the 
cleanest test of the Majorana nature of massive neutrinos indeed, Majorana mass terms 
violate lepton number by two units j^j giving rise to production or decay processes with 
AL = 2. Other mechanisms onderfymg the generation of nentnno masses, like the ones 
involving Higgs triplets |3| , can also provide a source of LNV. Currently, neutrinoless double 
beta (0u(3(3) nuclear decays (A, Z) — > (A, Z+2)e~e~ have become the most sensitive probe to 
search for the effects of very light Majorana neutrinos [4]. The underlying mechanism leading 
to these transitions is the conversion of two bounded neutrons in the initial nucleus into two 
bounded protons in the final one, making the knowledge of the nuclear wavefuntions the 
main limitation to achieve precise theoretical predictions. At the quark level, the elementary 
process dd — > uue~e~ is responsible for nuclear decays. 

The same simple mechanism would produce AL = 2 violation in hyperon decays, B^ — > 
B£l~l'~ , where L> 12 are hyperon states and 1,1' = e or /i. Examples of these decays are 
shown in Table 1: 
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TABLE I. Lepton number violating (AL = 2) decays of hyperons. The classification of these decays 
according to their change in strangeness (AS) is also indicated. 



Only one experimental upper limit of the channels listed in Table 1 has been reported 
so far, namely B(E~ — > p^T < 4.0 x 10 -8 l5j. A less restrictive AL = 2 decay mode 



in the charm sector has been reported in Ref. [6| with the following upper limit: L?(A+ — > 
< 7.4 x 10~ 4 . In the case of the decays listed in Table 1, two down-type (d 



2 



or s) quarks convert into two up-quarks changing the charge of hyperons according to the 
AQ = AL = +2 rule , as is shown in Figure 1. These quarks conversion are assumed to 
occur at the same space-time location and, therefore, they are driven by local four-quark 
operators. Therefore, the study of the relatively simpler case provided by Qv{3/3 hyperons 
decays may shed some light on the approximations used to evaluate the hadronic matrix 
elements relevant for similar nuclear decays. 

In the present paper we study the hadronic matrix elements of four-quark operators 
taken between initial and final hyperon states in the framework of the MIT bag model [7]. 
We use the effective low-energy Lagrangian proposed in Ref. which underlies AL = 2 
semileptonic transitions as the ones shown in Table 1. This method provides an evaluation 
of the hadronic matrix elements that does not use the approximation based on the insertion 
of intermediate states by factorizing the four-quark operators into two quark currents. The 
later approximation is commonly used in the evaluation of the hadronic matrix elements in 
neutrinoless double-beta decays of nuclei 4\ and hyperons il [lo|. 



II. EFFECTIVE LAGRANGIAN AND HADRONIC MATRIX ELEMENTS 

The most general form of the low-energy effective Lagrangian that is relevant for LNV 
semileptonic hyperon decays was given in Ref. [8| (the superscript c labels the charge 
conjugated spinor): 

Q2 

- C/3/3 = — — {ci(uTid)(uTjd) + c 2 [(uTid)(uTjs) + (uTis)(uTjd)} + c 3 (uTis)(uTjs)} 

A/3/3 

x {d 1 (eF k e c ) + d 2 (fiT kf i c ) + d 3 (eT kf i c + JiT k e c )} . (1) 

Here is a mass parameter corresponding to the physics scale for these processes and 
Q, di are dimensionless coefficients which represent the interaction strengths for the different 
channels. The dimensionless IYs are combinations of Dirac gamma matrices and depend on 
the physical mechanisms involved. The parameters and Lorentz structures involved in Eq. 
(CQ) depend on the specific underlying model and will contain some unknown parameters |8| . 
In the present paper, for the purpose of illustration, we will assume that only vector-axial 
structures are involved in fermionic bilinears althought it is not very difficult to consider 
other Lorentz structures. 
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At the lowest order in the interaction Lagrangian of Eq. (GQ), the 0v/3(3 decays of hyperons 
require the evaluation of the hadronic matrix elements of four-quark operators. The ampli- 
tude for any of the decays listed in Table 1, which we denote as B^{p) — > B~![ (p')l~ (pi)l~ (P2) 
with I — e or jj, (letters within brackets denote the four-momenta), is given by: 

m(b^ B+rr) = p^adjx^ir , (2) 

where 

= [«(pa)7i*( 1 - l5)luU c (pi) - (pi <r> p 2 )\ 

= 2g flv u(p 2 )(l + -y 5 )u c (pi) (3) 

is the properly antisymmetrized leptonic tensor, Cj and dj are the corresponding coefficients 
of the operators in the Lagrangian ([1]) and u c (pi) denotes the charge conjugated of spinor 
u(pi); note that L^ u becomes a symmetric tensor after using the charge conjugation property 



of the leptonic current 



The hadronic matrix element is: 

X*r* 2 = (Bt(p')\(u^(l - l5 )D)(u lu (l - j 5 )D')\B;(p)) , (4) 

where D, D' denote down-type quarks d or s. One way to compute the hadronic matrix 
element is to insert, between the product of quark bilinear operators, a set of intermediate 
baryonic states with the appropriate quantum numbers. Usually, one has to truncate the 
calculation by including only a few intermediate states which are supposed to be the domi- 
nant ones (for example, the S° and A hyperons in this case). This was done in Refs. {9!. [ic| 




u' 



FIG. 1. Feynman diagram describing the local interaction of Eq. ([T]) which converts two down-type 
quarks into two up quarks and two leptons, dd' —> uu'l~l~. 
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FIG. 2. Feynman diagram for the E —> pi I decay induced by the loop effect of a Majorana 
neutrino. 

in a model where the effects of virtual Majorana neutrinos is considered (see Figure |2] for 
the specific case of S~ — > pl~l~ decay). Next, one needs to use a set of form factors to 
parametrize the matrix elements of weak currents at each vertex; this procedure introduces 
a mo del- dependent input in the calculations. This approximation is good as long as only 
a few intermediate states and the low-energy behavior of the form factors give the domi- 
nant contribution. Note however that the loop integration becomes divergent and requires 
the introduction of an ad hoc regulator which can be identified with some average distance 
between quarks inside the hyperon 

In this paper, we use the MIT bag model of baryons to compute the matrix element given 
in 01]). Let us first note that, given the specific structure of quark currents in PJ, we can 
write it as follows: 

(B 2 (p')\(V - A) a (V> - A'UB^p)) = u(p') {Tip - T%] u(p) , (5) 

where u(p) are Dirac spinors describing the free hyperon states, and T^(P, q) are second- 
rank tensors that depends upon P = p + p' and q = p — p'. After using Gordon identities for 
the vector and axial currents, the most general form of the vertices can be parametrized as: 

= hig a p + ih 2 o- a i3 
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h D ^ 8 i ■ ^ 9 D u I • n 
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AM 2 a ^ p AM 2 ^ P 4M 2 " ^ AM 2 

'''VgYTs + ^e a ^P^ l5 , (6) 



2M ' 10 ' 2M 



5 



where M = (m + m')/2, with m(m') the mass of the initial (final) hyperon state. The 
coefficients hi are g 2 -dependent form factors which depend on the specific B\ — > E>2 tran- 
sition. Similarly, the axial vertex can be obtained by means of the following replacement: 
T^g = r^g(/ij — > gi) x 75. The form factors hi,gi in the vector and axial vertices have 
all a common dimension of energy. The contributions proportional to q/M are suppressed 
and the terms containing the Levi-Civita tensors do not contribute to the decay amplitude 
because of the symmetric leptonic tensor in Eq. (2). 



III. FORM FACTORS IN THE MIT BAG MODEL 

For definiteness, let us consider the specific example of the X~ — > pe~e~ transition; in 
this case only the operators with coefficients C2 and d\ in Eq. ([1]) give a contribution. In the 
framework of the MIT bag model, the vector and axial components of the hadronic matrix 
element in the £~ — > p transition, Eq. (4), can be written as follows: 

XT^ P {V) = u{p')T V a ,u{p) 

= (p\[Mty uu + M s a p uu ] + [Q d a y uu + QSTII^r) , (7) 



J T L J a/3 



{p\[Nty uu + N s a d ^ uu ] + \ptr uu + p s a v uu \\^) , (8) 



where we have defined (latin indices a, b, c, d denote flavor labels): 

M M-,ac = J d 3 x $ a ( x ) laijb ( x ^ . [^ C (x) 7/3 ^(X)], 



Qi d r c = J d 3 x[Mx)w 5 Mx)} ■ IMxh^M*)] • (9) 

In the above expresions ipi(x) denotes the wavefunction of quark with flavor i in the MIT 
bag model which is calculated according to Ref. [7} and it is reproduced in the Appendix. 
In the non-relativistic limit for baryons, the only non- vanishing matrix elements turn out 



to be the following: 



ds—^uu i j\ r sd— >uu] 



= 7^(wi_i,Wi_i,i?) 
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• (12) 



(p| E E I s ) + (pi E °^ r ' + E i s ) 
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(13) 



In the previous expressions r*, denote, respectively, the isospin and £/-spin rais- 
ing/lowering operators acting over the quark states in position i within the spin-flavor 
wavefunctions of |£~) and \p) (see Ref. 12J). Similarly, a^j refer to the fc-th component 



of the spin operator acting on the quark state in position j in the hyperon spin-flavor 
wavefunction. On the other hand, the functions 1Z, S, T and U introduced in Eqs. (10)-(13) 
arise from the integration over spatial coordinates of the quark wavefunctions in the MIT 
bag model; they depend upon the bag model parameters as shown in the expressions given 
in the Appendix. 

The matrix elements in Eq. (lOV ffTBl can be readily evaluated by using the quark model 
spin-flavor wavefunctions 12| of the E~ and p states. An explicit calculation yields: 
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and it follows (hereafter we omit the superscript labels of the hadronic matrix elements) 

X 00 (y) = 2K(w 1 - 1 ,w 1 - 1 ,R). (14) 
4 

X Qk (A) = --5 fej3 «S(wi_i,«;i_i,i2) (15) 
4 

X k0 (A) = --5 k<3 S(wi-i,wi-i,R) (16) 

X jk (V) = 2 ^T{wx-x,wi-x,R) [8j, 3 8 k , 3 - (1 - (J i)3 )(l - <S M y' + *[2(-l)' +1 + (-l) fc ]] 

+ ^ fe ZY(wi_i,wi_i,-R)} . (17) 

On the other hand, taking the non-relativistic limit of Eqs. (5) and (6) we get the 
following expressions for the non-vanishing hadronic matrix elements: 

X 00 (V) = h 1 + h 3 + h 5 = h 

X ok {A) = -{g 2 - g 5 + 2ih 13 )8 kj3 = f 2 8 k>3 , 

X k0 {A) = (g 2 + g 3 + 2ih 13 )5 ki3 = f 3 5 k)3 , 

Xj k (V) = -5 jk hi + (ih 2 - h 16 )(8j tl 8 k> 2 - 8 j)2 8 k>1 ) = f±8j k + / 5 (^,i4,2 - 8 j)2 5 k ^) , (18) 
or equivalently, from Eqs. (17) and (18): 

fx = 2TZ(wx-x,Wx-x,R) , 
4 

h = h = -^S(wx-x,wx-x,R) , 
2 

/ 4 = -T(wx-x,Wx-x,R) + 2U(wx-x,wx-x,R) , 

f 5 = 2iT(wx-x,Wx-x,R) , (19) 

for the effective form factors /j. Note that, in the non-relativistic limit, all these form factors 
should be evaluated at zero momentum transfer (q 2 = 0). Thus, we are not able to provide 
their momentum dependence; however, as in the case of the beta decays of hyperons, we 
may expect that these g-dependent effects would affect the decays rates by at most 10~20% 



given that they are a SU(3) symmetry breaking scale of order q/M 131 ] . 

The contraction of Lorentz indices in Eqs. (2)-(3) leads to the following simple form of 
the decay amplitude: 

M(5T -> pe-e~) = p^c 2 dx2u{p 2 ){l + ^)u c ( Pl ) ■ u{p') [A + B l5 ] u{p) , (20) 

A/3/3 

where A = Ahx + h 3 + h 5 = f\ — 3/4 and B = Agx + + ge- In the non-relativistic limit 
described above the numerical evaluations of Eqs. ffl9|) obtained in the framework of the 
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MIT bag model lead to A = 3.56 x 10 5 MeV 3 and B = (note the simplified expressions 
given in the Appendix for the functions 71, S, T and U). We have used the numerical 
values of the bag model parameters given in Ref. 7J]. 

From the above decay amplitude we get the following expression for the differential decay 
rate: 



-Y4 / 7 \ 2 

T F _ 

f 32vr 3 m 3 V A, 



dY Gi, ( cid 



F l c 2"l \ / 2 o 2 



(q 2 - 2mf)(u max - u 



min\ 



X 



{A 2 [(m^ + m p ) 2 -q 2 } + B 2 [(m^-m p ) 2 -q 2 }} , (21) 



where u max — u mm = (a/1 — 4m^/g 2 ) • A^fm 2 , m 2 , q 2 ) and X(x,y,x) denotes the triangle 
function. By using rjj- = 1.479 x 10~ 10 s jl4|, and integrating numerically the differential 



rate in the range 4m 2 < q 2 < (ms — m p ) 2 , we get the branching ratio for the di-electron 
channel: 



B h ^(Z~ -> pes-) = (^] ■ (4.65 x 10~ 13 MeV 2 ) . (22) 

If the take for weak coupling coefficients Cj,dj ~ 0(1) and A^ > 100 GeV, the branching 
fraction turns out to be extremely suppressed: £> bag (X~ — > pe~e~) < 10" 



-23 



Equivalently, we can define the following ratio [8j: 

= r(S- ne~u e ) ■ (23) 

Using the result given in ([22]), we get R^ ev {^~) ~ (^i/A^) 2 • (4.6 x 10~ 10 MeV 2 ). If we 
compare this number with the estimate given in Ref. 8j, which is based on pure dimensional 
arguments, the result of our present calculation appears to be smaller by six orders of mag- 
nitude. On the other hand, the corresponding result obtained by using the approximation 
based on the insertion of intermediate baryon states in the baryon-neutrino loop of Figure 



2 is i^XXE") = 7.2 x 10- 18 |(m ee )| 2 MeV~ 2 (from Table 1 in Ref. (lOj). Using the current 



limit |(m ee )| < 1 eV, we get R^J eu (Tr) < 7.2 x 10~ 30 or, equivalently, a branching fraction 
Z? loop (£~ — > pe~e~) < 7.3 x 10~ 33 . This upper limit on R l °°J ev (Yj~) is smaller than the one 
obtained in our bag model calculation by ten orders of magnitude (assuming C2 = d\ ~ 0(1) 



and App > 100 GeV). Note however, that the result of Ref. [lOj depends strongly on the 



cutoff used to regularize the integral over the virtual neutrino momentum in the loop. 
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IV. CONCLUSIONS 

In this paper we have used the MIT bag model [?| to evaluate the hadronic matrix el- 
ements of four-quark operators required to compute the rates of 0v(3f3 decays of hyperons. 
These four-quark operators appears in the most general AQ = 2 low-energy effective La- 
grangian [8] that describes the AL = 2 lepton number violation in hyperon transitions. This 
method avoids the use of the approximation based on the insertion of baryon intermediate 
states, which requires the knowledge of form factors for hyperon beta decays at very high 
energy scales. To our knowledge, this is a novel method for direct calculations of hadronic 
matrix elements of four-quark operators. 

As an specific example, we have considered the S~ — > pe~~e~ lepton number violating 
decay and have computed the non-vanishing form factors of the £~ — > p transition in 
the non-relativistic limit using the spin-flavor wavefunctions of the hyperon states. Using 
reasonable values for the order of magnitude of couplings and mass scales of the underlying 
New Physics, we have computed the upper limit on the branching ratio which turns out 
to be of order 10~ 23 for the S~ — > pe~e~ decay. Althought this result turns out to be 
ten orders of magnitude larger that the calculation based on a model where this decay is 



induced by a loop of baryons and light Majorana neutrinos [10], it is still very small to 
be accessible to sensitivities reached by current experiments. It shows, however, that the 
calculations based on models involving loops of virtual neutrinos and the insertion of virtual 
intermediate hyperon states may underestimate the true branching fractions. 
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Appendix 

In the MIT bag model the eigenf unctions of quarks confined within a baryon, which is 
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assumed to be a spherical bag of radius R, are given by [7| (note that k = ±1): 



^ nKm {x,t) = -= i/ 5 i™ \ e ~ lW 7 > ( 24 ) 



where w nK satisfies the eigenvalue condition tanu> nK = w nK /(l + w nK ), j|( K ±i) are spher- 
ical Bessel functions and U m are two- component Pauli spinors. These eigenfunctions are 
normalized according to 



/ 



(x,t) = 5 

nn' Okk' ^mm' > 

(25) 



and the normalization factors are defined from the following integrals: 



I 



2 j -2 / ' • »\ ' 2lO nK -|- K 



r drj (w nK r/R) 



o ' " ' AN 2 (w nK ) w nK + k 

1 2w nK + 3 
AN 2 (w nK ) w nK + K 



R rMrjl{w nK r/R) = — i + 3/t . (26 ) 



The general form of the integrals involving the product of four eigenfunctions required in 
our calculations of hadronic matrix elements, see Eq. (9), are 



Sx^^X^L^X.t^^X^L^X.t) = JJ^S R),(27) 

^1234 

d*x^\ i (x } t)4> L2 (x,t)ilj[(x } t)G k 4> Li (x } t) = —^5 mim2 al 4 S(w niKl ,w n3K3 ,R), (28) 

iVi ~ - 



1234 



d 3 x^ L (x,t)ajtlj L2 (x,t)tlj[ (x,t)a k ip L4 (x,t) = — — x [afaf T(w niK1 , w nsK3 , R) 

JV1234 

"t"^mim2 ^m3m4^(^iiiKi) ^»i3K3; -R)] (^9) 



where we have introduced the following notation: iV 12 34 = N(w niKl )N(w n2K2 )N(w n3K3 )N(w n4K4 ), 

-^j TliKilTli) C/c = Umi&kUmj and S nfi = 5ni«2 ( ^«3n4 ( ^KiK2 ^K3K4 • 
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By inserting the solutions given in Eq. (|24|) in the previous results, one gets: 
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(32) 
(33) 



A numerical evaluation of these radial integrals leads to the values of the the form factors 
in Eqs. (19). 
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